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ABSTRACT 

The ESA Gaia spacecraft has two Shack-Hartmann wavefront sensors (WFS) on its focal plane. They are 
required to refocus the telescope in-orbit due to launch settings and gravity release. They require bright stars to 
provide good signal to noise patterns. The ccntroiding precision achievable poses a limit on the minimum stellar 
brightness required and, ultimately, on the observing time required to reconstruct the wavefront. Maximum 
likelihood algorithms have been developed at the Gaia SOC. They provide optimum performance according to 
the Cramer-Rao lower bound. Detailed wavefront reconstruction procedures, dealing with partial telescope pupil 
sampling and partial microlens illumination have also been developed. In this work, a brief overview of the WFS 
and an in depth description of the centroiding and wavefront reconstruction algorithms is provided. 

Keywords: Astromctry, Gaia, wavefront sensor, Shack-Hartmann, centroid, maximum likelihood, wavefront 
reconstruction 

1. GAIA, THE NEED FOR REFOCUSING 

The ESA Gaia mission will provide astrometry of a billion objects in the Galaxy with unprecedent precision and 
accuracy. In addition, intermediate resolution spectra will be obtained for millions of sources. More details on 
the mission general goals can be found elsewhere 

The payload is composed of two twin off-axis three mirror anastigmatic telescopes (TMA) with a rectangular 
pupil of 1.45x0.5 m feeding a common focal plan^ (see also Kohley et al., this volume). In addition to the 
powered surfaces, several plane mirrors are required, two for the pupil plane beam combiner and two for a 
common periscope. Fig. [T] provides an overview of the overall optical system. Most of the payload, including the 
mirrors, focal plane and torus support structure are made of silicon carbide (SiC). This material combines low 
weight with high stiffness and thermal conductivity, providing a very homogenous temperature in-orbit. 
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Figure 1. Gaia payload overview (left) and unfolded telescope optical design (right). Courtesy Astrium. 



Figure 2. M2 Movement Mechanism (M2MM). Courtesy Astrium. 



Gaia will operate in the visible range (300-1050 nm) with a very high quality optical system (total wavefront 
error budget ~50 nm). The mechanical tolerances for such a folded TMA system are tight, and smaller that the 
typical perturbations estimated for the launch vibrations and gravity release. Focusing mechanisms have thus 
been incoroporated to move each secondary mirror (M2), the so called M2 Movement Mechanisms (M2MM), 
built by SENER. Each M2MM has a fully redundant set of actuators capable of orienting the M2 surface with 
five degrees of freedom (which is enough for a rotationally symmetric surface). Fig. [2] shows a model of the 
system. 

Detailed tolerancing simulations have been done at Astrium, with the M2MM as the sole compensator. They 
show that the design wavefront error can be recovered after actuation of the M2MM, with a very small residual 
contribution to the WFE (a few nm). 

2. THE GAIA WAVEFRONT SENSOR 

Two Shack-Hartmann WaveFront Sensors (WFS), built by TNO, are located on the Gaia focal plane to provide 
the information required to drive the M2MM. The structure of the WFS is made of invar, while the optical 
surfaces are made from fused silica. Both materials provide a good thermal match to the SiC CCD support 
structure. The optical design of the WFS is displayed in Fig. |3] It is based on an input slit (12"x30"), an 
spherical collimator, a microlens array (387 /im pitch, 378 /im diameter), a beamsplitter cube and two fold 
mirrors. Each Gaia output pupil is sampled with an array of 3x 10 microlcnses. The fold mirrors are required to 
project the WFS image plane at the location of the Gaia astrometric field image plane. That is, the WFS pattern 
is projected onto a CCD sharing the same focal plane than the other astrometric detectors (although mounted 
on reverse). The focal length of the microlenses is 18.5 mm, six times smaller than that for the collimator (111 
mm), so the WFS CCD is operated with a TDI period six times larger. The zero wavefront error reference is 
provided by a set of three optical fibres, which mimic perfect point sources in the entrance slit thanks to the 
beamsplitter cube. They are fed by a SLED diode with 825 nm central wavelength. Additional information on 
the WFS is available in the literature."^ 

3. CENTROIDS, CRAMER-RAO AND MAXIMUM LIKELIHOOD 

The Shack-Hartmann wavefront sensor provides information on the wavefront slopes measuring the centroid 
displacement of each lenslet image with respect to a given zero error reference. This process critically depends 
on the centroiding precision achieavable. The higher the relative precision, the fainter the observed star can be. 
There is a limit on the maximum precision achievable for a centroiding algorithm: the Cramer-Rao lower bound. 
For a PSF or LSF observed with CCDs, it can be expressed aPMSl 
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Figure 3. Gaia wavefront sensor. Top: schematic layout of the WFS. From left to right, the gaia exit pupils and WFS 
entrance slit, collimator, microlens array and focal plane can be seen. Middle: optical design layout. The second fold 
mirror (FM2) deflects the incoming light to a beamsplitter (BS), where the reference signal from the optical fibres can be 
injected. The first fold mirror (FMl) sends the light to the collimator Ml, which produces an image of the output pupil 
at the location of the microlens array (MLA). Finally, the lenslet images are projected onto the same focal plane used for 
the Gaia astrometric field. Bottom: mechanical overview of the WFS. Courtesy TNO.^^ 



where cr^ is the Cramer-Rao centroiding uncertainty on the rj axis of an image, Sk the discrete values of 
the PSF or LSF, in units of electrons, measured for each of the k samples of the image (CCDs can be binned 
on-chip, and each sample can be larger or equal than one physical pixel), S[. is the derivative of Sk with respect 
to the axis rj sample coordinate, r the read-out noise (in electrons) and b the homogeneous sky background (in 
electrons) . The units of CTj, are pixels. They can be converted to angles multiplying by the pixel size and dividing 
by the telescope focal length. 

This equation has extensively been used to assess the performance achievable by Gaia. It is worth noting that 
only the largest contribution is from those samples for which the PSF or LSF slope is higher. The information 
contained in samples far from the PSF quickly goes to zero. Real or simulated PSFs can be used to estimate the 
Cramer-Rao centroiding precision, The derivatives can be estimated by nearest neighbours finite differences: 



S'k - ; (2) 

Reaching the Cramer-Rao precision is not trivial. It has been shown that a good way of doing it is using 
maximum likelihood estimators.^ In this approach, the sample data {Sk} is fitted by a mathematical model 
Nk{xi, X2, ...jXri), which depends on several variables. One or two variables are centroid coordinates, for LSFs 
and PSFs, respectively. The other variables can be considered "shape" parameters (e.g. total flux, PSF width, 
microlens diameter, ...) required to provide a good agreement between the model and the observations. 

The problem of retrieving the centroid has thus been transformed in obtaining the set of parameters {xi, X2, ■■■■,x 
providing the best fit. Maximum likelihood methods are those minimising the quantity 
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where Sk — Nk is the difference between the actual value and the model for each sample k and y/r"^ + b + Sk 
is the measurement error for Sk- This is just a weighted least squares minimisation problem. Many non- linear 
numerical methods are available to carry out such optimisation (e.g. Gauss-Newton, Levenberg-Marquardt). 
Some of them have been implemented in public domain libraries. In this work, the Java Apache Commons- 
MatfQ implementation of the Levenberg-Marquardt algorithm has been used. 



4. MATHEMATICAL FORMULATION 

All least squares minimisation algorithms need an estimate of the derivatives of the PSF model with respect to 
each fitted variable: 



^ = ^ {xi,X2,-.-,Xn) (4) 

Some algorithms approximate the derivatives applying finite differences on a sufficiently small interval. They 
are better suited for complex PSFs difficult to model. On the other hand, purely analytical models with an 
explicit expression for the derivatives, can be much faster. The applicability of one approach or the other 
depends on the complexity of the problem. 

Well designed Shack-Hartmann sensors are typically diffraction limited. This means that the monochromatic 
PSF can be approximated by an Airy pattern. In the following, a collection of formulae will be given to estimate 
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the mimbcr of electrons collected by a sample N^,- and its derivatives with respect to five parametc;rs: the centroid 
coordinates x, y, the total flux and the lenslet diameter D and focal length /. Both the monochromatic and 
polychromatic cases will be considered. The value of these formulae, whose derivation is straightforward, is to 
have all of them together to facilitate the implementation by an interested reader. 

4.1 Monochromatic PSF for a single microlens 

The monochromatic irradiance Airy pattern I{x, y) produced by a diffraction limited lenslet is given by 

I{x, y, A) = lip) = lo (^^) ' = /o( Jo(p) + Mp)f (5) 
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where J„(r) is the n-th order Bessel function of the first kind, p is the reduced radial coordinate, xq and j/o 
the location of the PSF centre on the focal plane, D the diameter of a microlens, / the focal length of a microlens 
and A the wavelength. 

4.2 Monochromatic flux normaUsation: number of electrons collected 

The monochromatic PSF formula given above is only valid for small angles. However, it quickly goes to zero 
and can thus be integrated over the whole focal plane to obtain the total number of electrons N 
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where the following orthogonality relation for Bessel functions has been used 
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and the units of /q ai'e electrons per unit area 

^0 = ^ (9) 
4.3 Derivatives with respect to monochromatic PSF parameters 

The derivatives of the intensity I{x,y,X) with respect to the five variables xq, yo, D, f and Jq that define the 
location, width and intensity are given below. Properties of the Bessel functions have been used to remove p 
from the denominator and avoid singularities near the origin. 
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^{x,y,X) = -Io{Mp) + J2{p)){Ji{p) + Js{p))^ (12) 



gjix,y,X) = -loiMp) + J2{p)){Ji{p) + Up))y (13) 
— {x,y,X) = {Mp) + Mp)f (14) 

4.4 Non-SI units 

The following non-SI units have been used to make the interpretation of the results easier: The total number 

of electrons collected N instead of the peak intensity /g, the lengths are measured in samples (aJgampici J/sampic)! 
instead of metres {x, y). The derivatives with respect to these derived magnitudes can be computed using the 
chain rule. 
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4.5 Image composition and number of electrons per sample 

The monochromatic images including the contribution for all the Icnslets in the WFS will be modeled as the 
addition of npupii = 3 x 10 components per telescope and a diffuse light homogeneous background. Partially 
illuminated lenslets have not been considered, because they are difficult to model. 

npupii 1 

-'totai(a;,y,A) = ^ 7i(a;,j/, A) + sky (18) 

4.6 Number of electrons per sample: integration 

The monochromatic number of electrons N^-{i,j,\) collected at the sample (i,j) and wavelength A, and its 
derivatives ^^'"ag''^''^^ with respect to any variable ^, are computed integrating on the pixel rectangular bound- 
aries: 
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4.7 Polychromatic calculation 

The expressions above are valid for monochromatic or quasi-monochromatic hght sources, such as the SLEDs 
feeding the WFS caUbration optical fibres. However, an additional integration in wavelength is needed to produce 
a WFS stellar image. 
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where w{X) is the weighting function, F{X) is the stellar spectrum, in units of photons per unit area per unit 
time, T(A), is the Gaia + WFS spectral response funtion, R{X) the reflectivity of a single protected silver mirror, 
9 the number of mirrors in the telescope- WFS optical path (6 for the telescope and 3 for the WFS) and QE(A) 
the CCD quantum efficiency. The wavelength dependence of the MLA overall transmission has been neglected. 
The integration in wavelength has been approximated by a sum over discrete wavelengths 
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where the averages are carried out in nx intervals of equal length. 

5. SIMULATIONS: IMAGE GENERATION AND CENTROIDING PRECISION 

An image simulator and analyser has been developed at the Gaia SOC to produce artificial WFS images and 
test the fitting algorithm presented. It has been coded in Java, using existing libraries, such as those already 
developed by the Gaia Data Processing and Analysis Consortium (DPAC).'^ The least squares optimiser has 
been taken from Apache Commons Math, while the Bessel functions are computed with the Colt librarjj^ 

Pixel integrations use a single interval Gauss-Legendre algorithm with 3 AL x 4 AC knots per pixel sampling. 
Monochromatic and polychromatic calculations are possible. The former are appropriate for the calibration 
sources, while the latter can be used for stellar sources. The polychromatic calculations cover the wavelength 
range 300-1050 nm, which has been divided into 10 sub-intervals. A reference G2V spectrum has been used to 
simulate stellar sources. Poisson shot noise and Gaussian read-out noise can be added to the perfect noiseless 
images to simulate real observations. 

A set of 1000 monochromatic PSFs have been simulated with different centroid locations and noise. The 
total number of electrons is 10000, similar to the number of electrons needed to provide a good sampling of the 
gaia wavefront error for a perfect algorithm. The central SLED wavelength of 830nm has been used. They have 
subsequently been fitted by the maximum likelihood centroiding algorithm. The fit residuals are shown in Fig. |4] 
Because of the rectangular pixel size, they have been converted to fim to provide a common comparison scale. 

It can be appreciated that the PSF sampling only has a minor impact: the AL centroiding precision is only 
~16% better, despite the pixel size being three times smaller. The centroid uncertainties found agree very well 
(within 5%) with the estimates provided by the application of Eq. [T]to the observed PSFs. The values obtained, 
however, are worse (~47% AL, ~70% AC) than those obtained applying the following equation 

- '^^ =0.129Mni (26) 
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Figure 4. Centroid residuals. Single image Gaia WFS centroid fitting resituals are shown. The performance is only slightly 
better along scan, in spite of the pixel size being three times smaller. This shows that maximum likelihood centroiding 
can provide good results even with slightly undersampled images. The results agree withing 5% with the Cramer-Rao 
estimations carried out on the images according to Eq. [l] This demonstrates the method used is, indeed, optimal. 

which is the limit for a PSF sampled with infinitesimally small pixels and negligible read-out noiseEIIS It is 
worth noting that this expression provides a 30% better precision than that obtained for a gaussian function 
with a width equal to the Airy disc core. That is, significant centroiding information is stored in the Airy disc 
secondary maxima. 

Three conclusions can be drawn from this experiment. First, the centroiding method used is indeed optimal, 
and no further precison can be gained using more sophisticated algorithms. Second, simple rules of thumb for 
well sampled pure Poisson noise image models are not directly applicable to obtain the centroiding precision 
achievable under conditions of PSF undersampling or significcant read-out noise. Third, accurate predictions on 
the performance of a real wavefront sensor can be obtained computing the Cramer- Rao limit according to Eq. [l] 
on simulated images. 

In addition to simple PSF modelisation, full telescope pupil mono- and polychromatic simulations have been 
run to demonstrate the applicability to scenarios similar to real operations. The full WFS combined image is 
the sum of the Airy patterns corresponding to each lenslet. 

For a 3 X 10 sub-pupils image (partially illuminated lenslets are not considered), a total of 121 variables have 
to simultaneously be fitted by the algorithm. Two tests (monochromatic and polychromatic) have been carried 
out to assess if this is, in fact, possible. The lenslet images are arranged in a perfect 3 x 10 network, with a 39 
pix AL and 13 pix AC unit cell. The total signal is 10000 electrons per lenslet. Fig. [5]shows an illustrative WFS 
synthetic image. The polychromatic calculations have used the G2V stellar template. 

Both the mono- and polychromatic calculations converge. The time required for this brute-force optimisation 
is 8.8 min and 80 min, respectively for a regular desktop Core2 Duo at 3 GHz. However, the computation speed 
can significantly be increased applying the following recipes. 



Figure 5. WFS full pupil pattern. 3 x 10 sub-pupil polychromatic images are located in a perfect rectangular network, 
with a 39 pix AL and 13 pix AC unit cell. 

1. Fitting each lenslet image independently. This is a reasonable assumption for most wavefront sensors and, 
in particular, for Gaia, where the overlap between neighbour subimages is below 1 e~. It can be done 
either solving 3x10 different problems or limiting the extension of the PSF and its derivatives. 

2. Using Look Up Tables (LUTs) to store values for the PSF values and its derivatives. Interpolating inside 
LUTs is much faster than Bessel function evaluation and integration, and can be generalised to different 
PSF models. 

3. Using modified least squares algorithms where LUTs are also used to store the matrix used to estimate the 
parameter update for each iteration. 

All the suggestions above have been implemented for the Gaia WFS. The times required to fit a pattern 
have thus been reduced to around a minute for a single pattern applying suggestions 1 and 2 and down to a 
ms for a single image applying all the optinisations explained. It also seems that using the massive paraellel 
computing power of modern FPGAs or GPUs, maximum likelihood algorithms could be applied to the much 
more demanding task of real time closed loop adaptive optics wavefront sensing. 

6. WFS WAVEFRONT RECONSTRUCTION 

Typical commercial Shack-Hartmann wavefront sensors for on-ground applications have hundreds or thousands 
of lenslets. When these devices, providing a dense pupil sampling, are carefully aligned to the pupil and focal 
plane axes, the wavefront can easily be reconstructed. In fact, manufacturers usually provide ad-hoc custom 
software providing either the raw wavefront or its expansion into a polynomial series. 

The case for Gaia, however, is more complicated. On the one hand, only a small number of lenslets (3x10 
fully illuminated lenslets) per pupil is available. By design, the lenslet pitch was chosen to be a quarter of the 
pupil short dimension. This guarantees that 3x10 elements are always fully illuminated, although they will not 
cover the full pupil. This choice makes the device much less sensitive to microlens alignment errors or launch 
settings, at the price of leaving a significant portion of the telescope pupil unsampled. In addition, the Gaia 
WFS microlens array will only be rougly aligned with the telescope output pupil and the CCD. 

The objective for the Gaia in-fiight WFSs is to measure low order aberrations due to launch setting and 
gravity realease, so a sparse pupil sampling is enough. However, some usual simplifying assumptions are now 
not possible. In particular, the real location of each lenslet with respect to the telescope pupil and focal plane 
must be determined and considered. In addition, each lenslet subtends a significant portion of the output pupil, 
and they may not be considered infinitesimally small. 

A custom wavefront reconstruction algorithm has been developed. It takes many ideas and definitions from 
previous work by Astrium. The wavefront is expressed as a series of bidimensional Legendre polynomials, the 
first of which are given in Table [l] In addition, it directly converts the slopes into the polynomial expansion 
without reconstructing the wavefront as an intermediate step. The microlenses are not considered point-like. 
Pupil integration is carried out in a sort of ray tracing matrix computation. This method could be useful to 
other missions with few elements WFS. 



Table 1. Bidimensional Legendre polynomials. The WFS decomposes the WFE into a series of Legendre polynomials. 
The first nine terms are given below. The first three have no impact on image quality (piston, tilt x and tilt y). 
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7. WAVEFRONT LEGENDRE DECOMPOSITION 

Three reference systems must be defined. First, the telescope pupil reference system IZpupniOpupn, x,y), where 
Opupii is located in the centre of the WFS pupil, x and y are parallel to the exit pupil minor and major axes, respec- 
tively, and are normalised coordinates: x,y G [—1, 1]. Second, the lenslet reference system 7?.iensiet(Oiensiet) a^', y') 
has axes parallel to the directions upon which the lenslcts are distributed. The and y' directions are approxi- 
mately parallel to the WFS microlens array minor and major axes x, y, respectively. Third, the CCD reference 
system 7?.ccD(0iensiet) 2;", y"), where the centroids are estimated. The transformation relations between them 
are as follows 



The transformation between 7?.pupii and T^-iensiet is 
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where 6 is the angle between the lenslet array and the WFS pupil, and Dy are the WFS pupil dimen- 
sions and (Aa;, Ay) are the offsets between the centre of the microlens array and the WFS pupil in normalised 
coordinates. The transformation between TJ-icnsict and TZccu is 
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where 6' is the rotation angle between the microlens array and the CCD. 

For a given field of view, the wavefront error (WFE) is a bidimensional function of the WFS normalised pupil 
coordinates {x,y). 

WFE = WFE(a;, y) ^ ajLj {x, y) (29) 

j = l,2,...,ni (30) 

where the wavefront is expressed as a series of bidimensional Legendre polynomials, ni, = 9 is the number of 

2-D Legendre polynomials used. Legendre polynomials arc particularly convenient because they are orthogonal 
in the range [-1,1]. This means that the RMS WFE of the Legendre expansion is: 



RMS = V(WFE2) _ (WFE)^ 
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where dxj and dyj are the unidimensional Legendre orders in x and y and the sum is extended to all terms 
other than piston and tilt: d^j + dyj > 1. 

The WFS samples the wavefront surface at n locations in the exit pupil, corresponding to each lenslet. In 
fact, the WFS measures, for each lenslet i — l,2,...n, the centroid location differences {Sx'-,5y'/) between the 
images produced by the aberrated telescope and a reference zero WEE source (calibration optical fibre). Eor a 
lenslet z, the centroid displacement {Sxi,Syi) is the average of the wavefront error slope over the lenslet pupil 
area multiplied by the lenslet focal length /: 
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where the angular brackets in Eq. [34] denote average over a given subpupil of the the derivative of a given 2D 
Legendre polynomial. The number of points used in the average is equivalent to a number of rays "traced" for 
each polynomial. The lenslet integrals can reformulated as 
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and can be expressed in matrix form 
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where a is a vector with the {aj} coefficients of the Legendre expansion, that is the final output of the wave- 
front reconstruction process, that forms an overdetermined system of linear equations together with the matrix 
B and the centroid displacements vector 6. Solving that system of equations is equivalent to reconstructing the 
wavefront. B^ is the pseudoinverse matrix and can be computed using e.g. singular value decomposition (SVD) 
techniques. 

This method has the advantage that it could, in principle, be applied to partially illuminated lenslets just 
clipping the rays outside the telescope pupil from the average over the 2D Legendre polynomials. There are 
several ways of sampling an optical pupil with discrete rays. The most popular ones are the rectangular grid 
(default scheme in Code V) and the Gaussian quadrature polar grid (default scheme in Zemax). However, none 
of these methods provide an equal areas sampling scheme free of border effects for a circular pupil. The custom 
approach used in this work is described in the following. 

The unit circle is divided into annuli. Each of them is divided into equal area annular sectors. The radial 
and angular dimension of each sector is made as similar as possible. For each annular sector, the ray passing 
through the center of mass is used. For a given number of annuli riR in the unit circle, with approximately equal 



radial extension, the number of sectors for the annulus i that makes the radial and angular dimensions most 
similar can be estimated as: 

rii = round(27r«) (39) 
The area Ai for a sector in annulus i, with inner and outer radii r^-i, is 
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If the area for all sectors is equal, it can be shown that 
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Finally, the radial center of mass rem of an annular sector limited by ri_i and in an annulus with sectors 
is given by 
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Fig. [6] shows an example of circular pupil partition using five annuli. The centre of mass of the 94 annular 
sectors is shown on the unit circle. It is clear that the sampling is quite homogeneous. 



8. SIMULATIONS: WAVEFRONT RECONSTRUCTION 

Wavefront generation and fitting capabilities have been added to the Java package developed for Gaia WFS data 
analysis following the approach described before. The method has been tested in two different scenarios. 

1. The first one uses a wavefront composed of just the nine Legendre polynomial terms shown in Table [T] 

2. The second one is similar to the first one, but adding a very large RMS WFE, to simulate the (mostly) 
low frequency post-launch effects. 

For the first set of simulations 1000 different wavefront 2D Legendre expansions have been generated. Each 
one is the sum of the first six non-trivial Legendre terms shown in Table [ij Each term is scaled by a random 
gaussian factor. The total WFE has been set to the whole low frequency error budget of ~43.9 nm. Noiseless 
microlens image centroid displacements have been calculated and fitted using the algorithm presented before. 
The results obtained (see Table [2]) are very good, and the reconstruction error is negligible. 

The second simulations mimic the post-launch misalignment using the same approach above, but increasing 
the RMS WFE to a pessimistic value of 1 ^vci. The results are also shown in Table [2j It is clear that the Gaia 
WFS is capable of correctly identifying and measuring these aberrations. The telescopes could, in principle, be 
realigned using a single day of WFS observations. 
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Figure 6. Circular pupil partition. Five annuli have been used. The centre of mass of the 94 annular sectors covering the 
unit circle is shown. They provide a quite homogeneous sampling. Rays outside the telescope exit pupil are not clipped 
in the wavefront reconstruction process. This is the default sampling scheme used in this work. 



9. CONCLUSIONS 

A maximum likelihood algorithm has been developed to determine the centroid of the sub-pupil images in the Gaia 
Shack-Hartmann wavefront sensors. Each subpupil is assumed to produce a perfect Airy diffraction patterns. 
The formulation for both the mono- and polychromatic cases is presented. The centroid location precision 
achived reaches the maximum performance possible dictated by the Cramer-Rao lower bound. Full pupil 3 x 10 
images have been successfully fitted, demonstrating the applicability to real-case scenarios. Further optimization 
is possible separating the problem for each lenslets and using look-up tables bot for the PSF and intermediate 
matrices. Millisecond performance has been demonstrated for a single lenslet on a desktop computer. 

A wavefront reconstruction algorithm for wavefront sensors with few lenslets not sampling whole rectangular 
apertures have been developed. The location of each lenslet in the telescope pupil is considered. Ray tracing 
integration over each leanslet is carried out to account for the significant size of each element as compared to 
the telescope output pupil. The method has successfully been applied to the Gaia case, where it has been 
demonstrated that large post-launch low order aberrations can easily be identified and expressed as a series of 
bidimensional Legendre polynomials. Finally, the algorithms presented will be applied to real WFS data in the 



Table 2. Legendre fitting average residual errors for two different scenarios. 1000 different simulations have been run 
for each case. First, only low frequency terms compatible with the Gaia error budget. Second, large low frequency 
aberrations (up to 1 fim) simulating bad post-launch condition. The units are nanometers. No noise has been introduced 
in the simulations. It is apparent that the wavefront reconstruction works flawlessly for low frequency aberrations. 



Method 


Li 


is 




L7 


Ls 


L9 


1. Low frequency 

2. After launch 


6.25e-14 
1.44e-12 


5.47e-14 
1.23e-12 


2.14e-14 
4.86e-13 


3.29e-14 
7.55e-13 


7.99e-14 
1.78e-12 


6.33e-14 
1.42e-12 



near future. 
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